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Abstract 

A new class of cylindrically symmetric inhomogeneous cosmological 
models for perfect fluid distribution with electromagnetic field is obtained 
in the context of Lyra's geometry. We have obtained two types of solutions 
by considering the uniform as well as time dependent displacement field. 
The source of the magnetic field is due to an electric current produced 
along the z-axis.Only Fvi is a non- vanishing component of electromagnetic 
field tensor. To get the deterministic solution, it has been assumed that 
the expansion 9 in the model is proportional to the shear a. It has been 
found that the solutions are consistent with the recent observations of type 
la supernovae and the displacement vector /3(t) affects entropy. Physical 
and geometric aspects of the models are also discussed in presence and 
absence of magnetic field. 

PACS: 98.80.Jk, 98.80.-k 
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1 Introduction and Motivations 



The inhomogeneous cosmological models play a significant role in understanding 
some essential features of the universe such as the formation of galaxies during 
the early stages of evolution and process of homogenization. The early attempts 
at the construction of such models have been done by Tolman [T] and Bondi [2J 
who considered spherically symmetric models. Inhomogeneous plane-symmetric 
models were considered by Taub [3J d] and later by Tomimura [5] , Szekeres [B] , 
Collins and Szafron [7|, Szafron and Collins [3]. Senovilla !9J obtained a new 
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class of exact solutions of Einstein's equations without big bang singularity, rep- 
resenting a cylindrically symmetric, inhomogencous cosmological model filled 
with perfect fluid which is smooth and regular everywhere satisfying energy 
and causality conditions. Later, Ruiz and Senovilla [10) have examined a fairly 
large class of singularity free models through a comprehensive study of general 
cylindrically symmetric metric with separable function of r and t as metric co- 
efficients. Dadhich et al. [11] have established a link between the FRW model 
and the singularity free family by deducing the latter through a natural and 
simple in-homogenization and anisotropization of the former. Recently, Patel 
et al. [T2] have presented a general class of inhomogeneous cosmological models 
filled with non-thermalized perfect fluid assuming that the background space- 
time admits two space-like commuting Killing vectors and has separable metric 
coefficients. Singh, Mehta and Gupta [13] obtained inhomogeneous cosmolog- 
ical models of perfect fluid distribution with electro-magnetic field. Recently, 
Pradhan et al. |14j have investigated cylindrically-symmetric inhomogeneous 
cosmological models in various contexts. 

The occurrence of magnetic field on galactic scale is a well-established fact 
today, and its importance for a variety of astrophysical phenomena is generally 
acknowledged as pointed out by Zeldovich et al. [T5] • Also Harrison [TH] suggests 
that magnetic field could have a cosmological origin. As a natural consequences, 
we should include magnetic fields in the energy-momentum tensor of the early 
universe. The choice of anisotropic cosmological models in Einstein system of 
field equations leads to the cosmological models more general than Robertson- 
Walker model [17] . The presence of primordial magnetic field in the early stages 
of the evolution of the universe is discussed by many [H] — [27] • Strong magnetic 
field can be created due to adiabatic compression in clusters of galaxies. Large- 
scale magnetic field gives rise to anisotropies in the universe. The anisotropic 
pressure created by the magnetic fields dominates the evolution of the shear 
anisotropy and decays slowly as compared to the case when the pressure is held 
isotropic [28 , 29] . Such fields can be generated at the end of an inflationary epoch 
[5U] — [31] . Anisotropic magnetic field models have significant contribution in the 
evolution of galaxies and stellar objects. Bali and Ali [35] obtained a magnetized 
cylindrically symmetric universe with an electrically neutral perfect fluid as the 
source of matter. Pradhan et al. [36] have investigated magnetized cosmological 
models in various contexts. 

In 1917 Einstein introduced the cosmological constant into his field equa- 
tions of general relativity in order to obtain a static cosmological model since, 
as is well known, without the cosmological term his field equations admit only 
non-static solutions. After the discovery of the red-shift of galaxies and explana- 
tion thereof Einstein regretted for the introduction of the cosmological constant. 
Recently, there has been much interest in the cosmological term in context of 
quantum field theories, quantum gravity, super-gravity theories, Kaluza-Klein 
theories and the inflationary-universe scenario. Shortly after Einstein's general 
theory of relativity Weyl [37] suggested the first so-called unified field theory 
based on a generalization of Riemannian geometry. With its backdrop, it would 
seem more appropriate to call Weyl's theory a geometrized theory of gravitation 
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and electromagnetism (just as the general theory was a geometrized theory of 
gravitation only), instead a unified field theory. It is not clear as to what extent 
the two fields have been unified, even though they acquire (different) geomet- 
rical significance in the same geometry. The theory was never taken seriously 
inasmuchas it was based on the concept of non-integrability of length transfer; 
and, as pointed out by Einstein, this implies that spectral frequencies of atoms 
depend on their past histories and therefore have no absolute significance. Nev- 
ertheless, Wcyl's geometry provides an interesting example of non-Ricmannian 
connections, and recently Folland [35] has given a global formulation of Weyl 
manifolds clarifying considerably many of Weyl's basic ideas thereby. 

In 1951 Lyra [39j proposed a modification of Ricmannian geometry by in- 
troducing a gauge function into the structure-less manifold, as a result of which 
the cosmological constant arises naturally from the geometry. This bears a re- 
markable resemblance to Weyl's geometry. But in Lyra's geometry, unlike that 
of Weyl, the connection is metric preserving as in Remannian; in other words, 
length transfers are integrable. Lyra also introduced the notion of a gauge and 
in the "normal" gauge the curvature scalar in identical to that of Weyl. In con- 
secutive investigations Sen gD] , Sen and Dunn [UJ proposed a new scalar-tensor 
theory of gravitation and constructed an analog of the Einstein field equations 
based on Lyra's geometry. It is, thus, possible [3D] to construct a geometrized 
theory of gravitation and electromagnetism much along the lines of Weyl's "uni- 
fied" field theory, however, without the inconvenience of non-integrability length 
transfer. 

Halford [42] has pointed out that the constant vector displacement field 4>i in 
Lyra's geometry plays the role of cosmological constant A in the normal general 
relativistic treatment. It is shown by Halford 43J that the scalar-tensor treat- 
ment based on Lyra's geometry predicts the same effects within observational 
limits as the Einstein's theory. Several authors Sen and Vanstone [33] > Bhamra 
[43] , Karade and Borikar [35] , Kalyanshetti and Wagmode [37J , Reddy and Inna- 
iah [35J , Beesham [32] , Reddy and Venkateswarlu [50] , Soleng [5TJ , have studied 
cosmological models based on Lyra's manifold with a constant displacement 
field vector. However, this restriction of the displacement field to be constant 
is merely one for convenience and there is no a priori reason for it. Beesham 
[52] considered FRW models with time dependent displacement field. He has 
shown that by assuming the energy density of the universe to be equal to its 
critical value, the models have the k = — 1 geometry. Singh and Singh [53] — 
[55] . Singh and Desikan [57] have studied Bianchi-type I, III, Kantowaski-Sachs 
and a new class of cosmological models with time dependent displacement field 
and have made a comparative study of Robertson- Walker models with constant 
deceleration parameter in Einstein's theory with cosmological term and in the 
cosmological theory based on Lyra's geometry. Soleng [5TJ has pointed out that 
the cosmologies based on Lyra's manifold with constant gauge vector <j> will 
either include a creation field and be equal to Hoyle's creation field cosmol- 
ogy [56j— [60j or contain a special vacuum field, which together with the gauge 
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vector term, may be considered as a cosmological term. In the latter case the so- 
lutions are equal to the general relativistic cosmologies with a cosmological term. 

Recently, Pradhan et al. [5T], Casama et al. [55], Rahaman et al. [55] . 
Bali and Chandani [53] , Kumar and Singh [53], Singh [55] and Rao, Vinutha 
and Santhi [57] have studied cosmological models based on Lyra's geometry in 
various contexts. With these motivations, in this paper, we have obtained exact 
solutions of Einstein's field equations in cylindrically symmetric inhomogeneous 
space-time within the frame work of Lyra's geometry in the presence and ab- 
sence of magnetic field for uniform and time varying displacement vector. This 
paper is organized as follows. In Section 1 the motivation for the present work 
is discussed. The metric and the field equations are presented in Section 2, in 
Section 3 the solution of field equations, the Section 4 contains the solution of 
uniform displacement field (/3 = /3o, constant). The Section 5 deals with the so- 
lution with time varying displacement field ((3 = (3{t)). Subsections 5.1, 5.2 and 
5.3 describe the solutions of Empty Universe, Zeldovich Universe and Radiating 
Universe with the physical and geometric aspects of the models respectively. 
The solutions in absence of magnetic field are given in Section 6. Sections 7 and 
8 deal with the solutions for uniform and time dependent displacement field. 
Finally, in Section 9 discussion and concluding remarks are given. 



2 The Metric and Field Equations 

We consider the cylindrically symmetric metric in the form 
ds 2 = A 2 (dx 2 -dt 2 ) + B 2 dy 2 + C 2 dz 2 , 



(1) 



where A is the function of t alone and B and C are functions of x and t. The 
energy momentum tensor is taken as has the form 



T? 



(p + p)UiU J + pgf + Ef , 



(2) 



where p and p are, respectively, the energy density and pressure of the cosmic 
fluid, and Ui is the fluid four-velocity vector satisfying the condition 



U Ui 



-1, u l %i = 0. 



In Eq. Ef is the electromagnetic field given by Lichnerowicz 



E 



3 _ 



hih 1 UiU J + -gf 



hihP 



(3) 



(4) 



where p is the magnetic permeability and hi the magnetic flux vector defined 

by 
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where the dual electromagnetic field tensor *Fij is defined by Synge [BS] 



* F ij = —^-UjkiF u . (6) 

Here F^ is the electromagnetic field tensor and e^f-i is the Levi-Civita tensor 
density. 

The co-ordinates are considered to be co-moving so that u 1 = = u 2 = u 3 and 
it 4 = 4. If we consider that the current flows along the z-axis, then F12 is the 
only non- vanishing component of Fij . The Maxwell's equations 

F { ij;k] = 0, (7) 
= 4vrJ\ (8) 

■j 

require that F12 is the function of x-alone. We assume that the magnetic per- 
meability is the functions of x and t both. Here the semicolon represents a 
covariant differentiation 




The field equations (in gravitational units c = 1, G = 1), in normal gauge 
for Lyra's manifold, obtained by Sen [4] as 

13 3 

% - 2 9ljR + ~ 4 9ij ^ k< ^ = - 87rT yi (9) 

where <pi is the displacement field vector defined as 

& = (0,0,0,/8), (10) 

where (3 is either a constant or a function of t. The other symbols have their 
usual meaning as in Riemannian geometry. 



For the line-element (p}, the field Eq. 
following system of equations 



with Eqs. © and (10]) lead to the 
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BC 



B\a C14 Ai / Bi C\ 
~W ~C~ ~A\W ~c 



0, (15) 



where the subscript indices 1 and 4 in A, B, C and elsewhere denote ordinary 
differentiation with respect to x and t respectively. 



3 Solution of Field Equations 

Equations ([TT)) - ([1~5]) are five independent equations in seven unknowns A, B, C, 
p, p, f3 and F\2- For the complete determinacy of the system, we need two extra 
conditions which are narrated hereinafter. The research on exact solutions is 
based on some physically reasonable restrictions used to simplify the field equa- 
tions. 



To get determinate solution we assume that the expansion 9 in the model is 
proportional to the shear a. This condition leads to 

*-(§)"• < 16 > 

where n is a constant. The motive behind assuming this condition is explained 
with reference to Thorne [70], the observations of the velocity-red-shift rela- 
tion for extragalactic sources suggest that Hubble expansion of the universe is 
isotropic today within w 30 per cent [TTJ [75] • To put more precisely, red-shift 
studies place the limit 

^<0.3 

on the ratio of shear, c, to Hubble constant, H , in the neighbourhood of our 
Galaxy today. Collins et al. [73] have pointed out that for spatially homoge- 
neous metric, the normal congruence to the homogeneous expansion satisfies 
that the condition 5 is constant. 



From Eqs. (TT])-([ni), we have 

A44 A\ A4B4 A4C4 B44 B4C4 C\\ B\C\ 
~A~~A? + AB + ^4C B BC = ~~C BC~ 

and 

& n Fi2 _ C44 C11 B44 Bn 
fiB 2 ~ C~ ~C~ ~B W 

We also assume that 



= K (constant) (17) 
(18) 



B = f(x)g(t), 

C = f(x)k(t). (19) 
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Using Eqs. HU) and (Q2J) in JTSJ) and lead to 

fc£ _ (2n - 1) ff 4 
fe _ (2n + l) g ' 



(20) 



(n - 1} T- n — "7"*' (21) 

//n - A 2 = ^/ 2 - (22) 

Equation (|20[) leads to 

k = cg a , (23) 

where a = §^tj and c is the constant of integration. From Eqs. (f2"Tj) and (f2"3")l , 

we have 

.9 r 

where 

^ = na(q - 1) + a ^ = if 



n(a-l) + l ' n(l-a)-l 
Equation ([22]) leads to 

1 



f = expl-K(x + x o yj, (25) 
where xo is an integrating constant. Equation (f24l) leads to 

3=^ + ^')^), (26) 



where b = + 1)N and ci, C2 are integrating constants. Hence from (I23|) 
and (|26|) . we have 

fc = c(c 1 e b * + c 2 e- b *)(^ T ' . (27) 

Therefore we obtain 

B = exp (Ik{x + x ) 2 ^j ( Cl e bt + c 2 e~ ht ) ^ , (28) 

C = exp (^K(x + x ) 2 ^j c ( Cl e bt + c 2 e~ ht ) , (29) 

A = a(cie b * + c 2 e- bt ) : ^ 1 , (30) 
where a = 2a, C3 being a constant of integration. 

After using suitable transformation of the co-ordinates, the model ([T|) reduces 
to the form 

ds 2 = a 2 ( Cl e bT + cae-^^WT 1 ( d X 2 - dT 2 ) + e KX \c x e hT + c 2 e- bT ) VTU dY 2 
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+ e Kx2 ( Cl e bT + c 2 e- bT )WTTJdZ 2 , (31) 
where x + xq = X, t = T, y = Y , cz = Z . 

For the specification of displacement vector [3 within the framework of Lyra 
geometry and for realistic models of physical importance, we consider following 
two cases described in Sections 4 and 5. 



4 When (3 is a constant i.e. (3 = flo (constant) 

Using Eqs. ([28]). (129]) and ([30]) in Eqs. (fTTj) and (fT4|) the expressions for pressure 
p and density p for the model (13 lj) are given by 



8np = 



1 



2n(l-oQ 



K 2 X 2_ 2(3 + a)b 2 Cl c 2 



{I + l)^ 2 



where 



(2na 2 + a 2 + 2a - 2n + 3)b 2 ipf 



2(e + iy 



7 Pi 



8tt P 



2n(l-t») 



-3K 2 X 2 - 2K 



2b 2 (a - l)cic 2 



{2na 2 - a 2 - 2a - 2n + l)b 2 tpl 



2(1+1) 



(32) 



(33) 



yji = c ie bT - c 2 e~ bT , 

ijj 2 = Cl e bT + c 2 e- bT . 

From Eq. [T5)) the non- vanishing component F12 of the electromagnetic field 
tensor is obtained as 



jl b 2 {l-a) 



p2 

^12 - o_ //? , 1\2 



8tt (£+l) 2 



4(£+l)ciC2 + (l + a)V' 2 
4$ 



(34) 



From above equation it is observed that the electromagnetic field tensor in- 
creases with time. 



The reality conditions (Ellis [74] ) 

(i)p + p>0, (ii)p + 3p> 0, 

lead to 

& 2 (n-na 2 -l) V> 2 J^CiCa 2 

(i+i) 2 v| - (TTIM > 



(35) 
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and 



6 2 (4n - 4na 2 - a 2 - 2a - 5) tp 2 4b 2 (a + 5)cic 2 



[i + lf 



1$ {I + i)V^ 2 2 



'■'2 



respectively. 

The dominant energy conditions (Hawking and Ellis [75] ) 



(»)p-p>0, (n)p+p>0, 



lead to 



6 2 (a + l) 2 V 2 4& 2 (a + l)c lC2 3 



(^ + i) 2 v^I + 



I ' 2 



and 



fe 2 (n — na 2 — 1) ^ 2 4& 2 cic 2 



(^ + i) 2 v 2 2 + 



> K(KX 2 + 1), 



(36) 



+ o#>a> 2 >2if(2JO^ + l), (37) 



(38) 



respectively. The conditions ([56]) and ([3"T)l impose a restriction on constant 
displacement vector /3q. 



5 When j3 is a function of £ i.e. = 

In this case to find the explicit value of displacement field /?(£), we assume that 
the fluid obeys an equation of state of the form 



P = IP, 

where 7(0 < 7 < 1) is a constant. Using Eqs. 
(fTTj) and (fill) we obtain 



(39) 

([50")) and ([59]) in equations 



4^(1+7)^ : 



and 



2n(l-c) 



(l+7)/? 2 (t) 



, 9 46 2 cic 2 6 (n — na 2 — 1) ri 

-K 2 X 2 ~K 



(e + i)^ 2 (e + 1) 2 



K 2 X 2 (\ +7) + 2^7 



(40) 



3a 2 </> 2 



na — a 



2fr 2 cic 2 {(l-a)(l-7)-4} 

(e + i)4> 2 

1b 2 

-2a-2n+ 1)(1 + 7) - 2(na 2 - n+ l)}-y± 

r 2 



(41) 



Here we consider the three cases of physical interest in following Subsections 
5.1, 5.2 and 5.3. 
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5.1 Empty Universe 

Putting 7 = in (J35]) reduces to p = 0. Thus, from Eqs. (gDJ) and (gTJ), we 
obtain the expressions for physical parameters p and /3 2 (i) 



4np 



2n(l-c) 



Ab 2 c\C2 b 2 (n — na 2 — 1) ipf 



(£ + l)V 2 2 (£ + 1)2 VI 



(42) 



/3 2 (t) 



3a 2 ?/> 2 
From Eqs. (gSJ and 



<«+i) 



^a. 26 2 (a + 4) Cl c 2 6 2 (a + l) 2 ^ 2 



(£+l)V> 2 2 



(£+l) 2 ^ 



(43) 



we observe that p > and /3 2 (i) > according as 



b2{n - na2 -%l-K{KX 2 + l)rt> Ab2c ^ 



and 



(£+1)2 



K 2 X 2 ifj 2 



(l + l)' 



b 2 (a + l) 2 2 26 2 (a + 4) Cl c 2 



(44) 



(45) 



(£ + l) 2 rl ' (£ + 1) 
respectively. 

Halford [5] has pointed out that the displacement field <j>i in Lyra's geome- 
try plays the role of cosmological constant A in the normal general relativistic 
treatment. From Eq. (|43|) . it is observed that the displacement vector (3{t) 
is a decreasing function of time which is corroborated with Halford as well as 
with the recent observations [76J, [77] leading to the conclusion that A(t) is a 
decreasing function of t. 



5.2 Zeldovich Universe 

Putting 7 = 1 in Eq. ([39]) reduces to p = p. In this case the expressions for 
physical quantities are given by 



(3 2 (t) 



2n(l-c) 

3a 2 ^ 2 ~ Tmr 



8np = 87rp : 



K 2 X + K 



a 2 ^ 2 ~ Tmr 



4b 2 c lC2 b 2 (a + l) 2 t/jf 



(£+1)^ (^+1) 2 V 2 

46 2 cic 2 



-K 2 X 2 -K 



(l + l)^ 2 



b 2 (n - na 2 — 1) ipi 



(£+1)2 VI _ 

The reality condition (Ellis [71]) 

(i)/9 + p>0, (ii)p + 3p> 0, 

lead to 

b 2 {n~ na 2 — 1) 2 2 2 ^ 4fe 2 Cl c 2 



-Vr - K(ioc 2 + i)v 2 > 



(46) 



(47) 



(48) 
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5.3 Radiating Universe 

Putting 7 = | in Eq. ([3"9"| reduces to p = In this case the expressions for 
(3{t), p and p are obtained as 

2 



2r»(l-c) 



2K 2 X 2 + X 



26 2 (a + 5)cic 2 



3a 2 V> 2 

6 2 (na 2 - 2a 2 - 4a - n - 1) ^ 2 



(£+l)^ 2 



87175 



87rp 



2a 2 V> ; 



2r»(l- a ) 
+ 



2n(l-c) 

2aV 2 (£+1 ' L 



-if 2 X 2 - # 



3(* + l) 2 VI 

46 2 ciC2 6 2 (n — na 2 — 1) ?A 



(^ + i)^ 2 (i + iy ^ 



46 2 cic 2 b 2 (n — na 2 ~ \) ip\ 



(£+l)^ 2 



(£+l) 2 



(49) 
(50) 

(51) 



From Eq. ((49]), it is observed that displacement vector (3 is decreasing function 
of time and therefore it behaves as cosmological term A. 
The reality conditions (Ellis [74] ) 

(i)p+p>0, (ii)p + 3p> 0, 



are satisfied under condition 

The dominant energy conditions (Hawking and Ellis [75]) 



(i)p-p>0, (n)p + p>0, 



lead to 



b 2 (n — na 2 — 1) 

(7+iT 2 



^\ - K(KX 2 + 1)^| > 



4b 2 cic 2 



(52) 



Some Geometric Properties of the Model 

The expressions for the expansion 6, shear scalar cr 2 , deceleration parameter q 
and proper volume V 3 for the model ([31]) are given by 



b{n(l - a) + (1 + a)} i/>i 



(1 + l)aij^^ ^ ' 



a 1 



b 2 [{n(l - a) + (1 + a)} 2 - 3n(l - a)(l + a) - 3a] ^ 2 



2„2,/ , + ^2 



3(£ + l) 2 a 2 V 2 

6cic 2 (£ + l) 
n(l — a 2 ) (cie hT — c 2 e 



(53) 

(54) 
(55) 
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V 3 = V^9 = a 2 ^ 2 < e+1) e KX . (56) 
From Eqs. ([52]) and ||S1J we obtain 

a 2 {nil - a) + (1 + a)} 2 - 3n(l - a 2 ) - 3a 

7T7 = i — " ttt^-i r—r, — " = constant. (57) 

9 2 3{n(l - a) + (1 + a)} 2 v ; 

The rotation u> is identically zero. 

The rate of expansion Hi in the direction of x, y and z are given by 

A 4 nb(l — a) ipi 



Hy 



A (£ + 1) V2 



B 4 b V 



y B {£ + 1)^2 

Hz ~c (JTT)^- (58) 

Generally the model pip represents an expanding, shearing and non-rotating 
universe in which the flow vector is geodetic. The model (|3ip starts expanding at 
T > and goes on expanding indefinitely when ^t^jt < 0. Since j = constant, 
the model does not approach isotropy. As T increases the proper volume also 
increases. The physical quantities p and p decrease as F\ 2 increases. However, 
if ^fepp > 0, the process of contraction starts at T > and at T = 00 the 
expansion stops. The electromagnetic field tensor does not vanish when b 7^ 0, 
and a / 1. It is observed from Eq. (|55|) that g < when ci > and c 2 > 
which implies an accelerating model of the universe. Recent observations of type 
la supernovae [76[ 177] reveal that the present universe is in accelerating phase 
and deceleration parameter lies somewhere in the range — 1 < q < 0. It follows 
that our models of the universe are consistent with recent observations. Either 
when c\ = or c 2 — 0, the deceleration parameter q approaches the value (— 1) 
as in the case of dc-Sittcr universe. 



6 Solution in Absence of Magnetic Field 



In absence of magnetic field, the field Eq. (J9J) with Eqs. (|2|) and (|10|l for metric 
(p} read as 



1 

A 2 



B 44 
~B~ 



C44 


+2( 


'B 4 


CA 


~c 




,B ' 


c ) 


i 


(Al 


A 44 


C44 




\A 2 


A 


~ ~c 




(A\ 


A44 


B44 


'A 2 


\A 2 


A 


~ ~B~ 



B 4 C 4 






BC 


BC _ 


= 8 


C n \ 
C J 


= 8np + 


V 


B 


= 8np + 


V 



S;rp+ |/? 2 , (59) 



(60) 
(61) 
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1 


B\\ Cn 








B X C X 


BaCa 

1 — 


A 2 


B C 




KB 


4) 


BC 


BC 



^p-\(i\ (62) 

L>J! ^i^ '-Tr) (63) 



(64) 



B C A \B C 
Eqs. dSOj) and ([61]) lead to 

-B44 fin C44 Cn 

Eqs. (HU) and ([Ml) lead to 



Eqs. flUJ) and dHSD lead to 
which on integration gives 



^ - ^ = 0. (65) 



h a^- = 0, (66) 

9 r 



g = (c 4 i + c 5 ) PTU , (67) 
where C4 and C5 are constants of integration. Hence from (|23p and (|67[) . we have 

k= c(c 4 t + c 5 )T^> . (68) 

In this case (fT7| leads to 

1 



/ = exp ^-K(x + x ) 2 j. (69) 
Therefore, we have 

B = exp Q# (a; + x ) 2 ^j (c 4 t + c 5 ) <^> , (70) 

C = cxpQx(x + xo) 2 ^c(c4i + c 5 )FTT), (71) 

A = a(c 4 f + c 5 )Ti 1 f^r j (72) 
where a is already defined in previous section. 

After using suitable transformation of the co-ordinates, the metric (p} reduces 
to the form 

ds 2 = a 2 (c 4 T) 2 "^\dX 2 - dT 2 ) + e Kx2 (ciT)^!) dY 2 

+ e KX \ciT)T^TT)dZ 2 , (73) 
where x + xq = X, y = Y, cz = Z, t + ^ = T. 

For the specification of displacement field (3{t) within the framework of Lyra 
geometry and for realistic models of physical importance, we consider the fol- 
lowing two cases given in Sections 7 and 8. 
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7 When j3 is a constant i.e. j3 = (3q (constant) 



Using Eqs. (jT0"|) - (fT2")) in Eqs lCT)) and (f6"2)) the expressions for pressure p and 
density p for the model ([75]) are given by 



87175 : 



o 2 (c 4 T) 2 OTsr 



n(l — a 2 ) + a 1 1 



87175 : 



3n(l-a) 



a 2 (c 4 T)" 



(a + 1) 2 
n(l-a 2 ) + a) 1 



4 Po, 



^-K(2 + 3KX 2 ) 



(a + l) 2 

The dominant energy conditions (Hawking and Ellis [751) 

(*)p-p>0, (n)p+p>0, 



lead to 



and 



^/3 2 a 2 (c 4 T) 2 <i+°)" ) > + 2iOT 2 ), 



n(l — a ) + a I 1 
T 2 



(1 + a) 2 

respectively. 

The reality conditions (Ellis [74] ) 



> isT(l + iOT 2 



lead to 



and 



(i)p + p>0, (M)p + 3p>0, 
n(l-a 2 )+a\ 1 



(1 + a) 2 
2[n(l-Q 2 )+a] 1 



— >K(l + KX 2 ), 



(1 + a) 2 

The condition (f?6")) and (|75|) impose a restriction on (3q 



(74) 



(75) 



(76) 
(77) 



(78) 
(79) 



8 When (3 is a function of t 

In this case to find the explicit value of displacement field (3{t), we assume that 
the fluid obeys an equation of state given by (|59")l . Using Eqs. ([70)1 - (|?2l and 
(|39|) in Equations (|59| and (|62| we obtain expressions for p(t) and (3(t) given 

by 



?7r(l + j)p 



a 2 (c 4 T) d+°> 



n(l — a 2 ) + a 
(a + 1) 2 



T 2 



2K(i + irx 2 ) 



. (80) 
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(l+ 7 )/3 2 (i) = ^r^T 

3a 2 (c 4 T) 



n(l - a 2 ) + a | (1 - 7) 



(a + 1) 2 J T 2 

(81) 



+ 2X7 + A"X 2 (l + 37) 
It is observed that p > and /3 2 (t) > according as 

> K(l + KX 2 ), (82) 



n{l-a 2 )+a\ 1 . , ,. v . 



and 

'n(l-a 2 ) + al 1 



(a + 1) 2 J T 2 

> X 2 X 2 , (83) 



(a + 1) 2 J T 2 
respectively. 

It is worth mention here that by putting 7 = 0, 1, | in Eqs. (|5D|) and (fBTj) , 
one can derive the expressions for energy density p(t) and displacement vector 
(3{t) for empty universe, Zeldovich universe and radiating universe respectively. 
It is also observed that these three types of models have similar properties as we 
have already discussed above. Therefore, we have not mentioned the expressions 
for physical quantities of these models. 

Some Geometric Properties of the Model 

The expressions for the expansion 8, Hubble parameter H, shear scalar a 2 , 
deceleration parameter q and proper volume V 3 for the model (|73p in absence 
of magnetic field are given by 

0_ o „_ n(l- a ) + (l + q) 1 

a(l + a)c™ T a+=) 

2 _ {n(l-a) + (l + a)} 2 -3n(l-a 2 )-3a 1 

° n(l-a) 2n(l-c) + 2(l + a) V 00 ) 

3a 2 (l + a) 2 c 4 <1+Q) T (1+a) 

^-' Unil'TAl + aV (86) 

, , „ 2n(l-a)+(l+a) 

V" 3 = = a e (c 4 T) ~^ . (87) 



= constant. 



From Eqs. flU} and |(85j) we obtain 

cr 2 _ {n(l - a) + (1 + a)} 2 - 3n(l - a 2 ) - 3a 
¥ ~ 3{n(l-a) + (l + a)} 2 



The rotation to is identically zero. 

The rate of expansion Hi in the direction of x, y and z are given by 

A4 _ rt(l - a) 1 
A _ (1 + a) T' 
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H, 



y 



B 



(l + a)T' 



1 1 



C 



(l + a)T~ 



a 1 



(89) 



The model (|73p starts expanding with a big bang at T = and it stops expanding 
at T = oo. It should be noted that the universe exhibits initial singularity of the 
Point-type at T = 0. The space-time is well behaved in the range < T < Tq. 
In absence of magnetic field the model represents a shearing and non-rotating 
universe in which the flow vector is geodetic. At the initial moment T = 0, 
the parameters p, p, /3, 0, a 2 and H tend to infinity. So the universe starts 
from initial singularity with infinite energy density, infinite internal pressure, 
infinitely large gauge function, infinite rate of shear and expansion. Moreover, 
p, p, P, 9, <r 2 and H are monotonically decreasing toward a non-zero finite 
quantity for T in the range < T < T in absence of magnetic field. Since j 
= constant, the model does not approach isotropy. As T increases the proper 
volume also increases. It is observed that for all the three models i.e. for empty 
universe, Zeldovice universe and radiating universe, the displacement vector (3{t) 
is a decreasing function of time and therefore it behaves like cosmological term 
A. It is observed from Eq. (j86|) that q < when a < 2n+i which implies an 
accelerating model of the universe. When a = — 1, the deceleration parameter 
q approaches the value (—1) as in the case of de-Sitter universe. Thus, also in 
absence of magnetic field, our models of the universe are consistent with recent 
observations. 

9 Discussion and Concluding Remarks 

In this paper, we have obtained a new class of exact solutions of Einstein's field 
equations for cylindrically symmetric space-time with perfect fluid distribution 
within the framework of Lyra's geometry both in presence and absence of mag- 
netic field. The solutions are obtained using the functional separability of the 
metric coefficients. The source of the magnetic field is due to an electric current 
produced along the z-axis. F\% is the only non- vanishing component of elec- 
tromagnetic field tensor. The electromagnetic field tensor is given by equation 
(|34[) . p, remains undetermined as function of both x and t. The electromagnetic 
field tensor does not vanish if b ^ and a / 1. It is observed that in pres- 
ence of magnetic field, the rate of expansion of the universe is faster than that 
in absence of magnetic field. The idea of primordial magnetism is appealing 
because it can potentially explain all the large-scale fields seen in the universe 
today, specially those found in remote proto-galaxies. As a result, the litera- 
ture contains many studies examining the role and the implications of magnetic 
fields for cosmology. In presence of magnetic field the model pip represents 
an expanding, shearing and non-rotating universe in which the flow vector is 
geodetic. But in the absence of magnetic field the model ([70]) is found that in 
the universe all the matter and radiation are concentrated at the big bang epoch 
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and the cosmic expansion is driven by the big bang impulse. The universe has 
singular origin and it exhibits power-law expansion after the big bang impulse. 
The rate of expansion slows down and finally stops at T — > oo. In absence of 
magnetic field, the pressure, energy density and displacement field become zero 
whereas the spatial volume becomes infinitely large as T — > oo. 

It is possible to discuss entropy in our universe. In thermodynamics the 
expression for entropy is given by 

TdS = d( P V 3 ) + p(dV 3 ), (90) 

where V 3 = A 2 BC is the proper volume in our case. To solve the entropy 
problem of the standard model, it is necessary to treat dS > for at least a 
part of evolution of the universe. Hence Eq. (f9T)]) reduces to 

TdS = Pi + (p+p)^ + ^ + ^ >0. (91) 
The conservation equation T? • = for fTJ leads to 

A4 i?4 C4 \ 3 n a 3 n 2 I A4 B4 C4 



Therefore, Eqs. (|9Tj) and (|92j) lead to 



3 3 J [„^4 C4 



2 0fk + -f\2-j + ^ + -£)<Q. (93) 

which gives to j3 < 0. Thus, the displacement vector j3{t) affects entropy because 
for entropy dS > leads to (3(t) < 0. 

In spite of homogeneity at large scale our universe is inhomogeneous at small 
scale, so physical quantities being position-dependent are more natural in our 
observable universe if we do not go to super high scale. This result shows this 
kind of physical importance. It is observed that the displacement vector j3(t) 
coincides with the nature of the cosmological constant A which has been sup- 
ported by the work of several authors as discussed in the physical behaviour of 
the model in Sections 5 and 8. In recent time A-term has attracted theoreticians 
and observers for many a reason. The nontrivial role of the vacuum in the early 
universe generates a A-term that leads to inflationary phase. Observationally, 
this term provides an additional parameter to accommodate conflicting data 
on the values of the Hubble constant, the deceleration parameter, the density 
parameter and the age of the universe (for example, see Refs. [7B] and [79"]). 
Assuming that A owes its origin to vacuum interaction, as suggested partic- 
ularly by Sakharov [80J, it follows that it would, in general, be a function of 
space and time coordinates, rather than a strict constant. In a homogeneous 
universe A will be at most time dependent [81]. In the case of inhomogeneous 
universe this approach can generate A that varies both with space and time. In 
considering the nature of local massive objects, however, the space dependence 
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of A cannot be ignored. For details, reference may be made to Refs. [52], [53], 
[84]. In recent past there is an upsurge of interest in scalar fields in general 
relativity and alternative theories of gravitation in the context of inflationary 
cosmology [55] ISH] |57] . Therefore the study of cosmological models in Lyra's 
geometry may be relevant for inflationary models. Also the space dependence 
of the displacement field f3 is important for inhomogeneous models for the early 
stages of the evolution of the universe. In the present study we also find f3(t) as 
both space and time dependent which may be useful for a better understanding 
of the evolution of universe in cylindrically symmetric space-time within the 
framework of Lyra's geometry. There seems a good possibility of Lyra's geome- 
try to provide a theoretical foundation for relativistic gravitation, astrophysics 
and cosmology. However, the importance of Lyra's geometry for astrophysical 
bodies is still an open question. In fact, it needs a fair trial for experiment. 
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